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Abstract Matrix multivariate Pearson type II-Riesz distribution is defined 
and some of its properties are studied. In particular, the associated matrix 
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eigenvalues distributions are obtained. 
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1 Introduction 

When a new statistic theory is proposed, the statistician known well about the 
rigourously mathematical foundations of their discipline, however in order to reach 
a wider interdisciplinary public, some of the classical statistical techniques have 
been usually published without explaining the supporting abstract mathematical 
tools which governs the approach. For example, in the context of the distribution 
theory of random matrices, in the last 20 years, a number of more abstract and 
mathematical approaches have emerged for studying and generalising the usual ma¬ 
trix variate distributions. In particular, this needing have appeared recently in the 
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generalisation, by using abstract algebra, of some results of real random matrices to 
another supporting fields, such as complex, quaternion and octonion, see Ratnarajah 
et al. [26j, Ratnarajah et al. [27|], Edelman and Rao |l2l |. Forrester [l5j, among many 


others authors. Studying distribution theory by another algebras, beyond real, have 
led several generalisations of substantial understanding in the theoretical context, 
and we expect that it is more extensively applied when a an improvement of its 
unified potential can be explored in other contexts. Two main tendencies have been 
considered in literature, Jordan algebras and real norrned division__algebras. Some 
works dealing the first approach are due to Faraut and Koranyi 


Casalis, and Letac [3[, Hassairi and Lajmi [191 ] . Hassairi et al. 


141 1. Massam [241 ]. 
21], Kolodziejek 

231 ]. and the references therein, meanwhile, the second technique has been studied 


2C 


Gross and Richards 


, Dfaz-Garcfa 


Q, 


Dfaz-Garcfa and Gutierrez-Jaimez 


Dfaz-Garcfa 


, among many others. 

In the same manner, different generalisations of the multivariate statistical anal¬ 
ysis have been proposed recently. This generalised technique studies the effect of 
changing the usual matrix multivariate normal support by a general matrix mul¬ 
tivariate family of distributions, such as the elliptical contoured distributions (or 
simply, matrix multivariate elliptical distributions), see Fang and Zhang j3] and 
Gupta and Varga [3]. This family of distributions involves a number of known 
matrix multivariate distributions such as normal, Kotz type, Bessel, Pearson type 
II and VII, contaminated normal and power exponential, among many others. Two 
important properties of these distributions must be emphasise: i). Matrix multi¬ 
variate elliptical distributions provide more flexibility in the statistical modeling by 
including distributions with heavier or lighter tails and/or greater or lower degree 
of kurtosis than matrix multivariate normal distribution; and, ii). Most of the sta¬ 
tistical tests based on matrix multivariate normal distribution are invariant under 
the complete family of matrix multivariate elliptical distributions. 

Recently, a slight combination of these two theoretical generalisations have ap¬ 
peared in literature; namely, Jordan algebras has been led to the matrix multivariate 
Riesz distribution and its associated beta distribution. Dfaz-Garcfa [f| proved that 
the above mentioned distributions can be derived from a particular matrix multi¬ 
variate elliptical distribution, termed matrix multivariate Kotz-Riesz distribution. 
Similarly, matrix multivariate Riesz distribution is also of interest from the math¬ 
ematical point of view; in fact most of their basic properties under the structure 
theory of normal j-algebras and the theory of Vinberg algebras in place of Jordan 
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algebras have been studied by Ishi |22j and Boutouria and Hassiri Q, respectively. 

In this scenario, we can now propose a generalisation of the matrix multivariate 
beta, T and Pearson type II distributions based on a matrix multivariate Kotz-Riesz 
distribution. As usual in the normal case, extensions of beta, T and Pearson type II 
distributions involves two alternatives, the matricvariate and the matrix multivariate 
version;!], see Diaz-Garcfa fl, Diaz-Garcia jsS, Q , Diaz-Garcia and Gutierrez-Jairnez 
and Dfaz-Garcfa and Gutierrez-Sanchez m- 


This article derives the matrix multivariate beta and Pearson type II distribu¬ 
tions obtained from a matrix multivariate Kots-Riesz distribution and some of their 
basic properties are studied. Section [2] gives some basic concepts and the notation of 
abstract algebra, Jacobians and distribution theory. The nonsingular central matrix 
multivariate Pearson type II-Riesz distribution and the corresponding generalised 
matrix multivariate beta type I distribution are studied in Section [21 Finally, the 
joint densities of the singular values are derived in Section [21 


2 Preliminary results 


A detailed discussion of real normed division algebras can be found in Baez [l| and 
Neukirch et al. [251 ] . For your convenience, we shall introduce some notation, 


although in general, we adhere to standard notation forms. 

For our purposes: Let F be a held. An algebra 21 over F is a pair (21; m), where 
21 is a finite-dimensional vector space over F and multiplication m : 21 x 21 — >• A is 
an F-bilinear map; that is, for all A £ F, x, y, z € 21, 

m(x,\y + z) = A m(x; y) + m(x; z) 
m{\x + y, z) = \m(x; z) + m(y; z). 

Two algebras (21; m) and ((£; n) over F are said to be isomorphic if there is an 
invertible map ^ : 21 —>• (£ such that for all x, y £ 21, 

4>(m(x,y)) = n((f>(x), (j>{y))- 

By simplicity, we write m{x\ y) = xy for all x, y £ 21. 


1 The term matricvariate distribution was first introduced Dickey |ll| . but the expression matrix- 
variate distribution or matrix variate distribution or matrix multivariate distribution was later used 
to describe any distribution of a random matrix, see Gupta and Nagar [dJ] and Gupta and Varga 
□[ , and the references therein. When the density function of a random matrix is written including 
the trace operator then the matrix multivariate designation shall be used. 
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Let 21 be an algebra over F. Then 21 is said to be 

1. alternative if x(xy) = (xx)y and x(yy) = ( xy)y for all x,y € 21, 

2. associative if x(yz ) = ( xy)z for all x,y,z € 21, 

3. commutative if xy = yx for all x, y € 21, and 

4. unital if there is a 1 € 21 such that xl = x = lx for all x € 21. 

If 21 is unital, then the identity 1 is uniquely determined. 

An algebra 2t over F is said to be a division algebra if 21 is nonzero and xy = 
Oat => x = Oat or y = Oat for all x, y 6 2t. 

The term “division algebra”, comes from the following proposition, which shows 
that, in such an algebra, left and right division can be unambiguously performed. 

Let 2t be an algebra over F. Then 21 is a division algebra if, and only if, 21 is 
nonzero and for all a, b € 21, with b ^ Oat, the equations bx = a and yb = a have 
unique solutions x, y € 21. 

In the sequel we assume F = 3ft and consider classes of division algebras over 3ft 
or u real division algebras ” for short. 

We introduce the algebras of real numbers 3ft, complex numbers £, quaternions 
Sj and octonions 0. Then, if 21 is an alternative real division algebra, then 2t is 
isomorphic to 3ft, £, Sj or 0. 

Let 21 be a real division algebra with identity 1. Then 21 is said to be normed if 
there is an inner product (•, ■) on 21 such that 

(xy,xy) = (x,x)(y,y) for all x,y € 21. 

If 21 is a real normed division algebra , then 21 is isomorphic 3ft, £, Sj or D. 

There are exactly four normed division algebras: real numbers (3ft), complex 
numbers (£), quaternions (Sj) and octonions (0), see Baez l|. We take into account 
that should be taken into account, 3ft, £, Sj and 0 are the only normed division 
algebras; furthermore, they are the only alternative division algebras. 

Let 21 be a division algebra over the real numbers. Then 21 has dimension either 
1, 2, 4 or 8. Finally, observe that 

3ft is a real commutative associative normed division algebras, 

£ is a commutative associative normed division algebras, 

Sj is an associative normed division algebras, 

D is an alternative normed division algebras. 
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Let £n,m be the set of all n x m matrices of rank m < n over 21 with m distinct 
positive singular values, where 21 denotes a real finite-dimensional normed division 
algebra. Let 2l nXTn be the set of all n x m matrices over 21. The dimension of 
2I nxm Qver - g / 3 mri ' Let A £ 2l nxm , then A* = A 1 denotes the usual conjugate 
transpose. 

Table [0 sets out the equivalence between the same concepts in the four normed 
division algebras. 


Table 1: Notation 


Real 

Complex 

Quaternion 

Octonion 

Generic 

notation 

Semi-orthogonal 

Semi-unitary 

Semi-symplectic 

Semi-exceptional 

type 


Orthogonal 

Unitary 

Symplectic 

Exceptional 

type 

UP(m) 

Symmetric 

Hermitian 

Quaternion 

hermitian 

Octonion 

hermitian 

'- J m 


We denote by 6m the real vector space of all S £ 21 mxm suc h that S = S*. In 
addition, let be the cone of positive definite matrices S £ 2l mxm . Thus, 
consist of all matrices S = X*X, with X £ Dn,m\ then is an open subset of &m- 
Let Dm consisting of all D £ 2l mxm , D = diag(di,..., d m ). Let be 

the subgroup of all upper triangular matrices T £ 2l mxm such that ty = 0 for 
1 < i < j < m. Let Z £ £n,m, define the norm of Z as ||Z|| = \/tr Z*Z. 

For any matrix X £ 2l nxm , dX denotes the matrix of differentials ( dxij ). Finally, 
we define the measure or volume element (dX) when X £ 2t nxm ,©m, Dm or Vr n . n , 
see Dfaz-Garcfa and Gutierrez-Jaimez [jJ and Dfaz-Garcfa and Gutierrez-Jaimez £j]. 
If X £ 2I” xm then (dX) (the Lebesgue measure in 2l nxm ) denotes the exterior 
product of the fdmn functionally independent variables 

n m (3 

(dX) = f\ f\ dx^ where dx^ = dx[^\ 

i=lj=l k =1 

If S £ 6m (or S £ X^-(m) with tu > 0, i = 1,..., m) then (dS) (the Lebesgue 
measure in &m or in X^ r (rri)) denotes the exterior product of the exterior product 
of the m(m — l)/3/2 + m functionally independent variables, 

m m (3 

(dS) = f\ dsn /\ f\ ds\f. 

«=1 i>jk =1 
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Observe, that for the Lebesgue measure (dS) defined thus, it is required that S G tym, 
that is, S must be a non singular Hermitian matrix (Hermitian definite positive 
matrix). 

If A € then (dA) (the Legesgue measure in Dm) denotes the exterior product 
of the f3m functionally independent variables 

n fi 

(d\) = A A dx f ] - 

i =1 k= 1 


If Hi G vi, n then 

m n 

(H*dHi) = A A h i dh - 

i= 1 j=i+1 

where H = (H^H^)* = (hi,..., h m |h m+ i,..., h„)* G tdP(n). It can be proved that 
this differential form does not depend on the choice of the H 2 matrix. When n = 1; 
A, j defines the unit sphere in 2l m . This is, of course, an (m — l)/3- dimensional 
surface in 2l m . When n = m and denoting Hi by H, (HdH*) is termed the Haar 
measure on UP(m). 

The surface area or volume of the Stiefel manifold Vm,n is 


Vol(< n ) = [ (HidHj) 


yrn^mrip/2 

rAn/3/2] ’ 


( 1 ) 


where rf n [a] denotes the multivariate Gamma function for the space 6m and is 
defined as 


T^ 


a = 


' Ae<p£ 


etr{—A} | A| 


a—(m—l)/3/2—l 


(dA) 


— 1)0/4 


nr[a-(*-l)/9/2], 


Z=1 


and Re(a) > (m — l)/3/2. This can be obtained as a particular case of the generalised 
gamma function of weight n for the space 6m with k = (k\, & 2 ,..., k m ) G taking 
k = (0,0,... ,0) G and which for Re(a) > (m — l)/3/2 — km is defined by, see 


Gross and Richards 




16j and Faraut and Koranyi jl3 |. 

etr{-A}|A| a -( m - 1 ^/ 2 -A K (A)(dA) 


j Ae<p£ 

_ n m{m- 1)0/4 


^T[o + fei-(i-1)^/2] 


i =1 


= [«]K[a], 


( 2 ) 


(3) 


6 




where etr(-) = exp(tr(-)), | • | denotes the determinant, and for A £ &m 


771—1 


g*(A) = |A m | fcm Yl \M\ ki 


ki ki+1 


( 4 ) 


2=1 


with A p = (a rs ), r,s = 1, 2,... ,p, p = 1 , 2 , ..., m is termed the highest weight vector, 
see Gross and Richards jig], Faraut and Koranyi (Tdl ] and Hassairi and Lajmi jl9j ] : 
And, [a]« denotes the generalised Pochhammer symbol of weight k, defined as 


w i = 


J(a - (i - 1)0/2)*, 


i= 1 


m 


n m(m- 1)0/4, F[o + **-(*- 1)0/2] 

2—1 


rt[a] 


r&fal 


where Re(a) > (m — 1)0/2 — and 

(a)i = a(a + 1) ■ ■ ■ (a + i - 1), 
is the standard Pochhammer symbol. 

Additional, note that, if k = (p,... ,p), then q K ( A) = |A| P . In particular if p = 0, 
then q K { A) = 1. If r = (ti, t 2 , ■ ■ ■, t m ), h > t 2 > ■ ■ ■ > t m > 0, then q K+r (A) = 
q K (A)q T (A), and in particular if r = (p,p,...,p), then q K+r { A) = q K+p { A) = 
|A| p g K (A). Finally, for B € T fj{m) in such a manner that C = B*B € ©,«, 
9 k(B*AB) = g s (Ck(A), and ^(B*' 1 AB" 1 ) = (q K ( C))- 1 g K (A) = q_ K (C)q K (A), 
see Hassairi et al. 21], 


Finally, the following Jacobians involving the 0 parameter, reflects the gener¬ 
alised power of the algebraic technique; the can be seen as extensions of the full 
derived and unconnected results in the real, complex or quaternion cases, see Faraut 
and Koranyi [14] and Dfaz-Garcfa and Gutierrez-Jaimez [ 7 ]. These results are the 
base for several matrix and matric variate generalised analysis. 

Proposition 2.1. Let X and Y £ £n,m be matrices of functionally independent 
variables, and let Y = AXB + C, where A £ £n,n, B € £m,m and C £ are 
constant matrices. Then 


{dY) = | A* A | m d/ 2 1B* B | mn P/ 2 (dX). 


(5) 
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Proposition 2.2 (Singular Value Decomposition, SVD). Let X £ Cn,m be matrix 
of functionally independent variables, such that X = WiDV* with Wi £ Vm,n, 
V € 11 ^(m) and D = diag(di, • • • , d rn ) £ di > • • • > d m > 0. Then 

m m 

(dX) = 2- m vr e JJdf (n_m+1)_1 JJ(d 2 - d^(dD)(V*dV)(WJdWi), (6) 

i=l i<j 


where 



0, p = 1; 
— m, (3 = 2; 
—2m, P = 4; 
—4m, [3 = 8. 


Proposition 2.3. Let X 6 £n, m matrix of functionally independent variables, 
and write X = V i T. where Vi £ Vm,n and T € T f/(m) with positive diagonal 
elements. Define S = X*X £ Then 


(dX) = 2- m |S| /3(n - m+1)/2 ” 1 ((iS)(V^Vi). (7) 


Finally, to define the matrix multivariate Pearson type II-Riesz distribution we 
need to recall the following two definitions of Kotz-Riesz and Riesz distributions. 
From Dfaz-Garcfa [b]. 

Definition 2.1. Let S € /r € £n, m and k = (k\,k 2 , ■ ■ ■ , k rn ) £ Ui m . 

And let Y £ £n,m and U{ B) £ T u(n), such that B = 77(B)* £7 (B) is the Cholesky 
decomposition of B £ &m- Then it is said that Y has a Kotz-Riesz distribution of 
type I and its density function is 

pmnp/2+Y,r=i k iTt[np/ 2] 
irmnfj /2 [nP/2 , k] IXI n P/ 2 \@\ m P/ 2 

x etr {-ptr [S _1 (Y - p)*0~ 1 (Y - /*)] } 

x q K Y - /x)*0- 1 (Y - p)fY(S)- 1 ] (dY) (8) 

with Re(n/3/2) > (m — l)/3/2 — k m ', denoting this fact as 


Y ~ /C^ m (K,M,0,S) 
From Hassairi and Lajrni [lJ ] and Dfaz-Garcfa 




Definition 2.2. Let E £ and k = {k\, k 2 , ■ ■ ■, k m ) € 5R m . Then it is said that 
V has a Riesz distribution of type I if its density function is 


pam+YJILi k i 
Tm[o, «]|E| a g K (E) 

for V € and Re(a) > (m - 


etrj-ZJE- 1 V} | V| g K (V) (dV) (9) 

- l)/3/2 — km ; denoting this fact as V~ 7 Zm (a, k, E). 


3 Matrix multivariate Pearson type II-Riesz distribu¬ 
tion 


A detailed discussion of Riesz distribution may be found in Hassairi and Lajmi [IS] 
and Dfaz-Garcia [ 4 ]. In addition the Kotz-Riesz distribution is studied in detail 
in Diaz-Garcfa Q]. For your convenience, we adhere to standard notation stated in 
Diaz-Garcfa Q, Dfaz-Garcfa jo| . 

Theorem 3.1. Let = Si ~ Tl^' 1 (y/3/2, k, 1), k € 3ft and Re(V/3/2) > 

—k; independent of Y ~ JCR,^ m (r, 0, I n , I m ) 7 Re(n/3/2) > (m — l)/3/2 — t m . In 
addition, define R = S I_1 / 2 Y where S = S\ + \\Y\\ 2 . Then S ~ IZ^ 1 {{v + mn)f3/2 + 
Y^LiUik,!) independent o/R. Furthermore the density of R is 


r^[n/I/2]rf [(z/ + mn)f3/2 + k + YaL 1 U] 
nh mn / 2 Tm[nP/2, r]T^ [z//3/2 + k] 


(l-||R|| 2 )^ /2+fc “ 1 g T (R*R) (dR), 

( 10 ) 

where (1 — ||R|| 2 ) >0; which is termed the standardised matrix multivariate Pear¬ 
son type II-Riesz type distribution and is denoted as R k,T,l, 0 , I n , I m). 

Proof. From definition 12.11 and 12.21 the joint density of S 1 and Y is 
oc 4 V/2+fe_1 etr {-P (si + ||Y|| 2 )} q T (Y*Y) (dsi)(dY) 


where the constant of proportionality is 

pv(3/2+k prnnp/2+Y.ZiU r^n/3/2] 
T^P/2 + k\ Tr mn P/ 2 Tm [n/3/2 , r] 

Making the change of variable S = S\ — ||Y|| 2 and Y = S I ^ 2 R, by ((SJ) 

{ds 1 ) A (dY) = s pmn/2 {ds ) A (dR). 
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Now, observing that S = S\ — \ | Y| | 2 = S (l — ||R|| 2 ), the joint density of S and R 
is 

oc (1 - || R || 2 ) /3l ' /2+fe ~ 1 s pv/2+k-i etr {— q T (sR*R) (ds)(dR). 

Also, note that 


q T (sR*R) = q r ((s^I^R^^O = q T (sl m ) q T (R*R) = s^ ti q T (R*R). 


From where, the joint density of S and R is 
p(.v+mn)P/2+k+J2iLi *» 


rf [{v + mn)/3/2 + k + Yn=i *»] 


etr {-/3s} s K™*)/i/2+*+E£i U-i ( ds ) 


x 


r^[n^/2]rf [{v + mn)/3 /2 + A: + Yj?=i U] 


(1-IIRII 2 ) 


2\ v/3/2+k—l 


q T (R*R) (dR), 


vr /3mn/2r^ l [ n /5/2, r]r'j [zz/3/2 + fe] 

which shows that S IZi' 1 ((v + mn)/3/ 2 + 1), and is independent of R, 

where R has the density (1101) . □ 

Corollary 3.1. Let R ~ 'Pxx'^rnxni 1 '^ r ; 1, 0, I n , I?n) and define 


C = p- 1 / 2 U{e)*RU('Z) + fi 

where U(B) € Ty(n), such that B = £Y(B)*^(B) is the Cholesky decomposition of 
B G © G S G p > 0 constant and p G £n,m is a matrix of constants. 
Then the density of S is 

oc (1 - ptr S-^C - p)*Q-\C - p)) ,,fi/2+k - 1 

x 9T [/d(S)*” 1 (C-/i)*0^ 1 (C-/r)^(S)' 1 ] (dS) (11) 
where (l — p tr S~ 1 (C — /r)*© _1 (C — p)) > 0; with constant of proportionality 

ri[nP/2]r( [{y + mn)P/2 -k- ZT=i *»] p mn ^ 2 -^^ u 
7bW2 r £ [ n/ 3/2, —7-]rf [i^/3/2 - jfe]|£|0 n / 2 |0|0 m /2 

which is termed the matrix multivariate Pearson type II-Riesz distribution and is 
denoted as C ^ll^mxntb "A Pi /T ■ 

Proof. Observe that R = p 1 / 2 {C — /z)Z^(E)” 1 and 


(dR) = p mn P/ 2 \-X\-P n / 2 \®\-P m / 2 (dC). 


The desired result is obtained making this change of variable in (1101) . 


□ 
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Next we derive the corresponding matrix multivariate beta type I distribution. 

Theorem 3.2. Let R ~ t, p, 0, I n , S), and define B = R*R € tym, 

with n>m. Then the density of B is, 

oc iBifr-m+iWa-^i - ptr'£~ 1 B) uf}/2+k - 1 q T {B)(dB), (12) 


where 1 


ptrH X B > 0; and with constant of proportionality 


if [{y + mn)(3 /2 + k + Y1T= l *i] p^ mn / 2+ £)S= 1 

r^[n/3/2,r]rf [z//3/2 + fc]|£|"0/V(E) 


B is said to /iar>e a nonstandarised matrix multivariate beta-Riesz type I distribu¬ 
tion. 


Proof. The desired result follows from (flOl) . by applying ((2]) and then (JT|); and ob¬ 
serving that q T (U(f£‘)*~ 1 ~BU(f£‘)~ 1 ) = q- T (E)q T (B). □ 


In particular if S = I m in Theorem 13.21 we obtain: 


Corollary 3.2. Let R ~ P n Kfx m (^ k, r, 1,0, I n . I m ), and define B = R*R E 
with n>m. Then the density ofB is, 


rf [(v + mn)f3/2 + k + YT=i 
T m[n/3/2, r]r^ [nfi/2 + k] 


gi(n-m+l)/3/2-l^ 


pti-BY^ /2+k - 1 q T (B)(dB), 


(13) 

where 1 — ptrB > 0. B is said to have a matrix multivariate beta-Riesz type I 
distribution. 


Remark 3.1. Observe that alternatively to classical definitions of generalised ma- 
tricvariate beta function (for s ym metric cones), see Diaz-Garcia ja], Faraut and 
Koranyi 13] and Hassairi et al. 20], defined as 

Bm[a,K]b, t] 


f 

|g|6-(m-l)^/2-l ( 

0<CS<Clm 


= / 

|p|6— (m—1)/3/2—1 




_ Trn[a, t\ 

Tm[a -I- b, k + t] ’ 

where re = (h,k 2 , ■ ■ ■, k m ) € 5k m , r = (ti,t 2 , ■ ■ ■ ,t m ) € Re(a) > (m — l)fi/2 — k m 
and Re(6) > ( m — l)/3/2 — t m . From Corollary 13.21 and Dfaz-Garcfa and Gutierrez- 
Sanchez [id . Theorem 3.3.1], we have the following alternative definition: 
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Definition 3.1. The matrix multivariate beta function is defined an denoted as: 


B*J[a,k-b,r] = / ^b-{ m -i)p/ 2 -i^ _ tr B) a+fc " V(B)(dB) 

J 1—tr B>0 

= / |F| fo -( m - 1 )/3/2- 1 (i _|_ fr F) _ ( a+m6+fc+ Stei*i)g r (F)(dF) 

•'Re<p£ l 

T^[a + k]Tm[b, t] 
rf [a + mb + k + YT=i **] 

Also, observe that, when m = 1, then r = t and k = k and 


[a, k ; b, f] 


rf[g + fc]r^[6 + t] 

rf [a + 5 + A: + f] 


i3 [a, fc; 6, t] 


Finally observe that if in results in this section are defined k = 0 and r = 
(0,... ,0), the results in Dfaz-Garcfa and Gutierrez-Jaimez [§| are obtained as par¬ 
ticular cases. 


4 Singular value densities 


In this section, the joint densities of the singular values of random matrix R ~ 
P n K fxm(b k > T > !. °, In, I m) are derived. In addition, and as a direct consequence, 
the joint density of the eigenvalues of matrix multivariate beta-Riesz type I distri¬ 
bution is obtained for real normed division algebras. 


Theorem 4.1. Let 5\,... ,5 m , 1 > <5i > • • • > 5 m > 0, be the singular values of the 
random matrix R P ii^nxm(b T > !> °, In, Im)- Then its joint density is 

m \ v/3/2+k-l 


r)m /3m 2 /2 +q ” 

^ ' TT ( e2 \ 

rt,[l3m/2}B*J[vP/2, k- nfi/ 2, r] * 


2 \ ( 71 — 771 + 1)^/2 — 1/2 


X 


n (j . 2 -$ 


\ i= 1 

/3 Cr (D 2 ) 


A dSi 


(14) 

i \ 7 ' / 

i<j V A nV \i=l / 

/or 1 — pYaLi $i > 0- WTiene £ is defined in Lemma UTiH D = diag(di,..., <5 m ), and 
Ck(-) denotes the zonal spherical functions or spherical polynomials, see Gross and 
Richards \l$] and Faraut and Koranyi \l\. Chapter XI, Section 3]. 


Proof. This follows immediately from (jlOlh First using ((6|), then applying H]) and 
observing that, from [itl Equation 4.8(2) and Definition 5.3] and Faraut and Koranyi 
[fill . Chapter XI, Section 3], we have that for L € fpm, 

c?( Z) = Cf ( lm) [ 9«(HZH*)(dH), 

JHeU^(m) 
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□ 


Finally, observe that Si = y / eig, i (R*R), where eigj(A), i = 1,... ,m, denotes the 
i-th eigenvalue of A. Let Aj = eigj(R*R) = eigj(B), observing that, for example, 
Si = y/Xi- Then 

m mm 

f\d6i = 2- m l[\; 1/2 

i= 1 i=l i=l 

the corresponding joint densities of Ai,..., A m , 1 > Ai > • • ■ > A m > 0 is obtained 
from (11411 as 

m / m \ v/3/2+k—l 


'll 


/3m 2 /2+g 


Tm[f3m/2]J3m l3 [v[3/2,k]nl3/2,T\ f =1 




i —1 


X 


n < A < - a 


Ct{G) 


i<j \ 

for 1 - YllLi A i > o, where G = diag(Ai,..., A m ) 


fr A". 

L m) \j = i 


5 Conclusions 


As visual examples, different Pearson type II-Riesz densities for m = 1 are showed 
in figures Q] and [2J 

Recall that in octonionic case, from the practical point of view, we most keep in 
mind the fact from Baez jdj, there is still no proof that the octonions are useful for 
understanding the real world. We can only hope that eventually this question will be 
settled on one wav or another. In addition, as is established in Faraut and Koranyi 
[14;] and Sawyer 28] the result obtained in this article are valid for the algebra of 
Albert , that is when hermitian matrices (S) or hermitian product of matrices (X*X) 
are 3x3 octonionic matrices. 
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Pearson ll-Riesz type density functions 



X 


Figure 1: With v = 15, n = 18 and t = 7 


Pearson ll-Riesz type density functions 



X 


Figure 2: With v = 3, n = 18 and k = 0 
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